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We constrain the possibility of a new pseudoscalar coupling between the muon and proton using 
a recent measurement of the 2S hyperfine splitting in muonic hydrogen. 


Recent measurements of 2 S — 2 P transition frequen¬ 
cies in the exotic atom constituted by a proton orbited 
by a muon HE] find the proton charge radius to be 7a 
smaller than the 2010-CODATA 3] value obtained using 
ordinary hydrogen and e— p scattering. The 2S hyperfine 
splitting deduced from the same measurements shows ex¬ 
cellent agreement with predictions [Tj . The discrepancy 
in the proton radius has generated a lot of interest, in¬ 
cluding the invocation of new fundamental interactions 
as an explanation. Here, we focus on the implications 
of the hyperfine splitting for new interactions between 
the muon and proton. Specifically, we consider the pos¬ 
sibility of a new pseudoscalar particle that couples to the 
muon and proton. Such an interaction is spin and ve¬ 
locity dependent and has a negligible effect on the Lamb 
shift (which is used to extract the proton radius) in the 
nonrelativistic limit j¥j , but has a significant effect on the 
hyperfine splitting. 

The measured value of the 2S hyperfine splitting 
(HFS) m 

AEjjfs = 22.8089 ± 0.0051 meV, (1) 

is to be compared with the theoretical prediction [5] 

A E% FS = (22.9843±0.0030) —(0.1621±0.0010)rz+i5£)j, 

( 2 ) 

in meV, where the Zenrach radius [B] 

r z = 1.045 ± 0.004 fm, (3) 

is obtained from e —p scattering^] SE a is the contribution 
to HFS from the new pseudoscalar interaction. Taking 
the experimental and theoretical uncertainties in quadra¬ 
ture, the best-fit to the experimentally measured A Ehfs 
and r z occurs for r z = 1.045 fm and SE a = —0.006 meV, 
and 

- 0.018 meV < SE a < 0.006 meV at 2a. (4) 

We now compute SE a , and subject it to the above 2 a 
constraint. In the nonrelativistic (NR) limit, the pseu¬ 
doscalar vertex becomes 

h = w(p0*75«(p) w^^iwX - *x ,f ^Tx, 


1 The use of the value of rz obtained from e — p scattering is 
appropriate here because the correction to rz from using the 
new ii — p interaction arises at loop order. 
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FIG. 1. Pseudoscalar exchange in muonic hydrogen. 


where x an d x' are 2-component Pauli spinors. The p—p 
interaction in terms of the muon line (given by Xyt c fu) 
and the proton line (given by x P > °>) is then (see Fig.nl), 

• (P - P')Xy TikfXpVp • ( P - P ')Xp , 

and the NR scattering amplitude for p + P —>• p' + P f 
is 


iM = , with q=p-p' = P' - p. 


The couplings of the light pseudoscalar a of mass m a to 
the muon and to the proton are f^ and f p , respectively. 
Then, 

M = - z ^r p x'j^-qxu x'J<7p-<lXp q4^ 
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with the relative angle averaged for the s wave. The 
effective Hamiltonian is 
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so that 
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FIG. 2. The values of f p f p allowed at 2er lie between the solid 
curves. 


where if is the wave function of the 2 S state: 


ip(r) 


1 

2 V /27r “B 


(1 



Here, as = —— is the Bohr radius for muonic hydrogen 
with m r = mpm p /{mp + m p ), the reduced mass of the 
system. On convolving, we obtain 

SB. = MAsAfp (iM , (5) 

Sm^nip o7t \m r J 

where 


F(x) = 1 — x 2 


a 2 * + 2x 2 
2 (a + a;) 4 


( 6 ) 


It is important to distinguish between m r and m p in the 
equations above. The m r dependence comes from the 
Bohr radius as, and m M from the NR reduction. The 
function F(x) interpolates between 1 and 0 for x = 0 and 
x —» oo which is consistent with decoupling behavior. 

In Fig. [2j we show the 2c allowed values of fpf p as a 
function of m a . The region between the solid curves is 
allowed. We restrict m a < 100 MeV so that f p f p remains 
comfortably in the perturbative regime. 

Note that in the potential model of the proton with 
nonrelativistic quarks, the proton pseudoscalar coupling 
f p arises from the pseudoscalar couplings f u , fd of the up 
and down quarks, which are of the same order of magni¬ 
tude. In this simplified picture, we have f P = |/d — \ f u 
(as for the magnetic moments). If f u = fd, we have the 
simple result, f p = f u = fd- 

In principle, the anomalous magnetic moment of the 
muon places a stringent independent constraint on f p 


since for the m a of interest, pseudoscalar couplings yield 
a negative contribution to a p [7], ;8] 0 while the measured 
value is higher than the standard model expectation: 
Aap = a e p xp - afp = (29 ± 9) X 10~ 10 0. However, 
the scalar sector may be more intricate than envisioned 
here, and may offer a fine-tuned (and perhaps unnatural) 
cancellation of the pseudoscalar contribution. 

For the sake of comparison, the QED contribution at 
leading order is 


d-HQED 


9r_ § 3 {r) 

6 4 m^rrip 


cr 


p • 


Here, (« 2), and g v (« 5.5857) are the gyromagnetic 
ratios for the muon and proton. Correspondingly, 
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The above QED result, though simple, represents the 
first three significant digits of the dedicated theoretical 
calculation, and is consistent with the recent measure¬ 
ment of Ref. II]. 

The ratio of the pseudoscalar contribution to the lead¬ 
ing QED contribution is 


SE a _ 2 fpfp ^ / m a \ 

SEqed 4:Tragpg p \m r ) 

2 fpfp ^ _ m 2 a 2 + 2 (m a /m r ) 2 ' 

47ra g p g p \_ m% 2(a + m a /m r ) 4 

In sum, the 2S hyperfine splitting in muonic hydro¬ 
gen constrains the product of the pseudoscalar cou¬ 
plings of the muon and proton fpf p to lie in the 
2(7 ranges [-0.00040,0.00013], [-0.00173,0.00058] and 
[—0.015,0.005] for m a = 0, 10 MeV and 100 MeV, respec¬ 
tively. As the pseudoscalar mass is further increased, the 
constraint is weakened. The couplings have no impact 
on the discrepant measurements of the proton radius. 

For m a < 100 MeV, no direct limits on f p f p exist from 
colliders, although limits for higher m a were obtained by 
the CMS experiment using the dimuon channel in pp col¬ 
lisions m- The CMS upper limits on the cross section 
times branching fraction, a ■ B{pp — > a — > p + p~), di¬ 
rectly constrain f p f p in the mass ranges, 5.5 — 8.8 GeV 
and 11.5 — 14 GeV. CLEO’s nonobservation of the decay 
J/il) —l ya, with a invisible pT] gives the 90% C.L. con¬ 
straint |/p| < 0.029 (assuming the J/ij) — a coupling to 
be f p ) for m a < 100 MeV [!J. 
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2 In Eq. (11) of Ref. [J], Cp should be replaced by |Cp| 2 , since 

Cp, as defined in Eq. (9) therein, is complex. 
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